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ABSTRACT 

We report here the first polarimetric measurements of the pulsars in the J0737-3039 binary neutron star system 
using the Green Bank Telescope. We conclude both that the primary star (pulsar A) has a wide hollow cone of 
emission, which is an expected characteristic of the relatively open magnetosphere given its short spin period, 
and that A has a small angle between its spin and magnetic dipole axes, 4 ± 3 degrees. This near alignment of 
axes suggests that A's wind pressure on pulsar B's magnetosphere will depend on orbital phase. This variable 
pressure is one mechanism for the variation of flux and profile shape of pulsar B with respect to the orbital phase 
that has been reported. The response of pulsar B to the A wind pressure will also depend on the particular side 
of its magnetosphere facing the wind at the spin phase when B is visible. This is a second possible mechanism 
for variability. We suggest that pulsar B may have its spin axis aligned with the orbital angular momentum 
owing to A's wind torque that contributes to its spindown. Monitoring the pulsars while geodetic precession 
changes spin orientations will provide essential evidence to test detailed theoretical models. We determine the 
Rotation Measures of the two stars to be -1 12.3 ± 1 .5 and -1 18 ± 12 rad m"^. 
Subject headings: pulsars: polarization - radiation mechanism: non-thermal 



1. INTRODUCTION 

The double pulsar system, PSR J0737-3039, recently re- 
ported by Burgay et al. (2003) and Lyne et al. (2004) is likely 
to unlock many mysteries concerning isolated neutron star 
magnetospheres and winds as well as the nature of the pul- 
sar emission mechanism. A 22.7-ms pulsar (A) and a 2.77 sec 
pulsar (B) revolve about each other in a 2.4h, nearly edge-on 
orbit. The two stars show interesting and contrasting emis- 
sion properties. Pulsar A exhibits a complex double profile 
structure. Apart from an eclipse when A is behind B (with 
respect to our sight line), this pulsar does not show significant 
variation in its flux. The flux and profile structure of B show 
remarkable variations as a function of orbital phase (Lyne et 
al. 2004). 

A simple pressure balance calculation leads to the conclu- 
sion that the MHD wind of pulsar A will surpress a large frac- 
tion of the quasi-static B magnetosphere. The eclipse dura- 
tion of pulsar A is roughly consistent with this idea (Lyne et 
al. 2004; Kaspi et al. 2004). Detailed calculations of the 
A-wind/B-magnetosphere interactions are required for com- 
parisons with observations (Arons et al. 2004). In this model 
a bow shock, magnetosheath, magnetopause and magnetotail 
structures, which are relativistic analogs of solar wind/Earth 
magnetosphere structures, are established around pulsar B. 

In this work, we first present polarimetric measurements on 
these two stars. These observations were conducted at the 
National Radio Astronomy Observatory^ Green Bank Tele- 
scope (GBT). After describing our observational setup in |2] 
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we present our results in gS] Our observations place strong 
constraints on the geometrical orientation of A's rotation and 
magnetic axes. We also present our Rotation Measure (RM) 
measurements for the two stars. In l^we discuss implications 
of our polarimetric results within the context of the Arons et 
al. model. 

2. OBSERVATIONS 

Observations were performed at the GBT with two back- 
ends: BCPM (Berkeley-Caltech Pulsar Machine; see e.g., 
Camilo et al. 2002) and GBPP (Green Bank-Berkeley Pul- 
sar Processor; see e.g.. Backer, Wong & Valanju 2000). The 
pulsar system was observed at four different frequencies on 
2003 December 11, 19, 23, & 24, and 2004 January 1. The 
bandwidths of our BCPM observations were 96 MHz at 2.2 
GHz and 1.4 GHz, and 48 MHz at 0.82 GHz and 0.43 GHz. 
The GBPP bandwidths were 28 MHz at 2.2 GHz, 1.4 GHz 
and 0.82 GHz and 1 1 MHz at 0.43 GHz. 

On 2003 December 1 1 we observed A with the GBPP at 
1.4 GHz using the published ephemeris (Burgay et al. 2003) 
and an integration time of 100s. The GBPP provides syn- 
chronously averaged, coherently dedispersed profiles that al- 
low calculation of all Stokes parameters. 

We first detected and then generated an ephemeris for 
B with the BCPM filter bank total-power data of 2003 De- 
cember 1 1 . This ephemeris was used for GBPP observations 
of B at 0.82 GHz on 2003 December 19. GBPP observa- 
tions continued with A at 0.43 GHz (briefly), 0.82 GHz and 
2.2 GHz and B at 0.43 GHz. We also measured a correlated 
(fully linearly polarized) noise source and observed pulsars 
with known polarization characteristics (PSR B0656-I-14 and 
PSR B0919H-06) for calibration purposes. GBPP observations 
of PSR B1929H-10 at 0.82 GHz by I. Stairs for another pro- 
gram on 2003 December 26 were also analyzed for calibra- 
tion. 

3. RESULTS 
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Fig. 1. — Polarization cliaracteristics of J()737-3039A at 0.82 GHz. The FiG. 2. — contours for tlie RVM fit for A . Tlie filled contouns repre- 

inset gives the two portions of the profile 'miiTor-folded' to show similarity sent Icr and 3a confidence Umits. The Une contours show the emission cone 

of outer edges. See text for details. opening half-angle p for each solution. 



The top panel of Fig. shows the average profile of A at 
0.82 GHz. The pulsar shows two prominent regions of emis- 
sion with sharp outer edges, and shallower inner edges. The 
dotted and dashed lines indicate observed average profiles in 
linear and circular polarization, respectively. We also show 
the linear polarization position angle as a function of pulse 
longitude in the bottom panel of Fig.^ The longitude region 
-100° to -70°, and probably the region from -1-90° to -1-105°, 
are affected by the presence of quasi-orthogonal polarization 
modes. 

The two emission regions display rough mirror symmetry 
in both total intensity and linear polarization. We show in the 
inset on the figure two portions of the profile "mirror-folded" 
to align the outer edges and emphasize this point. These mir- 
ror symmetries suggest to us that these two portions of the 
pulse arise from traversals of a wide-angle hollow cone at the 
same magnetic pole rather than the two opposite poles. We 
have chosen the arrangement of the profile in rotational lon- 
gitude with steep outer edges in Fig. Rowing to what is often 
seen in more compact hollow cones (e.g., Manchester 1996). 

Assuming the rotating vector model (RVM; Radhakxishnan 
& Cooke 1969), we performed a least squares fit to the posi- 
tion angle sweep to obtain the value of the angle (a) between 
/2 and SI and the angle (/3) between the line of sight and jl at 
closest approach. For this purpose, we ignored the longitudes 
mentioned above where the PA is likely corrupted by quasi- 
orthogonal emission modes. The fit was done via a simple 
grid search in a and /? (e.g. Nice et al. 2001). This allowed 
us to map out contours for the solutions, which are plot- 
ted in Fig. 121 Two solutions appear, one with a ^ 5° and one 
with a ~ 90°. The large-a solution is unrealistic since it re- 
quires a beam opening angle 2p ^ 180°, a geometry where 
the RVM is likely to be invalid. Furthermore, in this case 
we would expect to see emission from both magnetic poles, 
which is inconsistent with the mirror symmetric profile. We 
favor the small-a solution since it avoids both these problems 
(see Fig. |3lfor a diagram of this geometry). Unfortunately, 
this leaves (3 completely unconstrained. Measuring f3 would 
give a lower limit on the misalignment between fl and the or- 
bital angular momentum, since the orbital inclination is well 
known. However, in the smaU-a solution /3 ^ p, so we might 



expect l3 ~ 50° based on opening angles of other millisecond 
pulsars. 

The characteristics of B's emission are much more com- 
plex. Lyne et al. (2004) show that the amplitude and pulse 
profile vary with orbital phase. Our sensitive data exhibit lo- 
cal maxima in four different "windows" in the orbit (/ to IV^). 
Within the limit of our sensitivity, this 4-window behavior and 
details of profile evolution within the windows are frequency 
independent over 0.43 GHz to 1.4 GHz. The signal is seen 
in almost all other orbital phases at very low flux levels ex- 
cept for a phase range of 30° to 92°, where no flux is detected 
above 0. 1 % of that in window /. The centroid of the strongest 
emission (windows / & //) is ^ 30° prior to when A is behind 
B, and the centroid of the phase range when B is invisible 
is ^ 30° prior to when A is in front of B . Due to the pulse 
structure evolution, we present the results separately for the 
four windows in Fig. @]from our 0.82-GHz observation with 
the same conventions as in Fig. [2 In window 1 the pulse is 
polarized at the 15% level and shows orthogonal polarization 
modes, while in window 2 the polarization is less than a few 
per cent. A detailed study of individual pulses from B and 
their evolution as a function of orbital phase by Ramachan- 
dran et al. (2004) concludes that the subpulse fluctuations are 
similar to those found in other pulsars of similar period and 
pulse morphology. 

3.1. Rotation Measures 

We have measured the RM of the two pulsars from our 0.82 
GHz observations. Assuming linearity of polarization posi- 
tion angle sweep within the small bandwidth of 28 MHz, we 
have derived RM values for pulsar Aand pulsar B from win- 
dow /. We have obtained consistent values of -112.3 ±1.5 
rad m"^ and -1 18 ± 12 rad m"- towards A & B, respectively. 
These values are not corrected for ionospheric contributions. 
Within our measurement errors, we do not see significant or- 
bital phase dependent variations. 

4. DISCUSSION 

* with the convention followed by Lyne et al. (2004), the ranges of orbital 
phases (true anomaly) we have taken for the four windows are: 186° —216°, 
252° -300°, 336° -30° and 90° - 132°, respectively 



Axis Orientations of J0737-3039 Pulsars 
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Fig. 3. — The proposed geometry of the emission beam of PSR J0737- 
3039A, as viewed from above the pole. The grayscale circles indicate the 
emission cone, centered around the magnetic axis (jl). The rotation axis is 
given by f2. The 'dash' line indicates the line of sight trajectory around f2. 
The radiating vectors from /J are projections of field lines diverging from the 
pole which nominally give the PA of linear polarization. The two emission 
regions in the profile are marked as a and b. 



A remarkable feature of both the Lyne et al. (2004) and 
our GBT observations of J0737-3039B is the dependence of 
flux and pulse profile, including polarization properties, on 
orbital phase, but not on radio frequency. We concur in gen- 
eral with Lyne et al. that the variable influence of pulsar A on 
the magnetosphere of B is the likely source of the changes. 
In a companion paper (Arons et al. 2004) we develop a de- 
tailed model in which A's MHD wind confines the B magne- 
tosphere. We outUne this model below, and make coimections 
between observational features that both stimulated the model 
making and provide directions for future work. 

In our model the MHD wind from A decelerates through a 
relativistic bow shock that envelops B's magnetosphere. Be- 
tween this shock and the boundary of B's closed magneto- 
sphere (B's magnetopause) the shocked A wind plasma forms 
a relativistic hot layer whose optical depth to synchrotron ab- 
sorption at 500 MHz is at least 100, and may be as high as 
5000. This layer (the magnetosheath) is the likely cause of 
the ecUpse of A by B. In this model we further interpret the 
faintest region of B emission near its superior conjunction as 
the result of absorption by the same magnetosheath layer that 
in these phases is between the observer and the pulsar emis- 
sion region deep within the B magnetosphere. In Kaspi et 
al. (2004) we report that the A ecUpse is asynmietric with 
slower flux decrease on ingress and a more rapid recovery on 
egress. As discussed above the strongest and faintest regions 
of B emission are also asymmetric: both precede times of 
conjunction by ~ 30°. In our model these are attributed to 
the prograde rotation of B that leads to an asymmetric mag- 
netopause and magnetotail. The frequency independence of 
the eclipse light curve over our observing bands is attributed 
to sharp boundaries, high optical depths and partial covering 
in space and/or time. 

We favor the interpretation of our polarization observations 
of pulsar A that has the dipole field axis nearly aligned with 
the spin axis, which itself is oblique to the orbital plane and 
the line of sight. Pulsar B would then experience a signifi- 
cant difference in wind pressure and content around the orbit. 



This variable pressure will affect the size of the polar cap and 
current structure that could lead to the B variations. 

The A wind, asymmetric or not, will produce a propeller 
torque on B owing to B's rotation. This torque contributes 
to the observed P of B. In this model the B emission results 
from voltage and e±-avalanche current deep within its mag- 
netosphere similar to that in normal pulsars. The "normalcy" 
of B emission is supported by these polarization observations 
as well as the single pulse study described elsewhere (Ra- 
machandran et al. 2004). 

When we observe pulsar B, it presents a different face of its 
closed and confined magnetosphere to the A wind depending 
on the orbital phase. This variable orientation will also change 
the B magnetosphere and, we suggest, its polar cap and the 
beamed emission we detect. At this point we cannot distin- 
guish between these two possible sources of modification of 
the secondary star - variable A wind pressure impinging on 
B, and variable B internal structure at the rotation phase of 
observation. They both are seemingly consistent with the fre- 
quency independence of the reported phenomena. Both cases 
will lead to variations in the polar cap size and relativistic cur- 
rent structure of B, which will affect the observed flux and its 
pulse morphology. 

In the case of the variable wind pressure, the phase of mod- 
ulation of B is expected to precess around the orbit as the spin 
of A undergoes geodetic precession. Geodetic precession of 
B will have a more complex effect on its internal response to 
a wind pressure at the spin phase of observation owing to the 
multiplicity of angles involved. 

Our proposed geometry of A's beam will be testable as 
geodetic precession moves the observer through the cone and 
into a region of invisibility during the 70-year cycle (Lyne 
et al. 2004). Is it improbable that we are seeing both pul- 
sars at the same time? We suggest that the A wind torque 
on B has aUgned pulsar B's spin axis with the orbital angular 
momentum over time. In this case B would need to have its 
dipole axis at a 90°, and therefore we will continue to view 
B independent of geodetic precession. An interesting question 
is whether the torque of the A wind can misalign the dipole 
while aligning the spin. The full history of the evolution of 
the system and the dueling magnetospheric winds remains to 
be written. In its infancy B might have had a short spin pe- 
riod and strong magnetic field. If so, then the early wind from 
B could have dominated the A magnetosphere and altered its 
spin and magnetic dipole. 

4.1. Summary 

Our polarimetric observations indidate that pulsar A's spin 
and magnetic axes are nearly aligned. This leads us to con- 
sider two possible mechanisms for the variability of B with 
orbital phase - (1) a pole to equator asymmetry of the A wind, 
and (2) a rotational asymmetry of the force balance radius of 
the B magnetosphere as it adjusts to the A wind. The interac- 
tion also contributes to the spin-down torque on the B pulsar 
which may also lead to alignment of the B spin and orbital 
momentum vectors. Geodetic precession will provide an im- 
portant means of separating the relative importance of the var- 
ious effects in the coming years. Theoretical calculations are 
underway to explore these ideas quantitatively. 

We thank the GBT staff, and in particular Carl Bignell, 
Frank Ghigo, Glen Langston and Karen O'Neil, for extensive 
help with the observations and very useful discussions. We 
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